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Abstract 

The time dependent Schrodinger equation for 7r+-7r~ pairs, which are emitted from 
the interaction zone in relativistic nuclear collisions, is solved using wavepacket states. It 
is shown that the Coulomb enhancement in the momentum correlation function of such 
pairs is smaller than obtained in earlier calculations based on Coulomb distorted plane 
waves. These results suggest that the experimentally observed positive correlation signal 
cannot be caused by the Coulomb interaction between pions emitted from the interaction 
zone. But other processes which involve long-lived resonances and the related extended 
source dimensions could provide a possible explanation for the observed signal. 
PACS: 25.70.-z 



1 Introduction 



Originally, the method of intensity-interferometry ( "Hanbury-Brown and Twiss interferome- 
try", HBT) was developed in astronomy and applied to determine the angular radii of stars, 
which were too small to be resolved in ordinary optical instruments Shortly afterwards it 
was recognized that a similar procedure could also be applied to pions emitted in high energy 
reactions, whereby the size of the pion-production volume becomes a measurable quantity 
1^]. However, it is clear that the different physical environment in high energy- or heavy ion 
reactions asks for certain modifications of the original formalism. The following facts have 
to be taken into account: 

1. The wavelength A of visible light is about 500 nm and the coherence length of a photon 
A ~ r c with a typical emission time of r ~ 10~^ s is about 7 orders of magnitude larger. 
Therefore, the photon can be well described as a plane wave. A typical pion, produced 
via A(1232)-decay with a momentum po = 150 MeV/c has the de Broglie wavelength 
Ao ~ 8 fm, and the coherence length is defined as (|^, p. 158) A = A^/AA ^ 6 fm. 
Here, AA is the uncertainty in the wavelength, related to the energy width /S.E ~ 120 
MeV of the A(1232) resonance via AA = 2ti fi/lS.E. Since Ao and A are of the same 
size, the pion cannot be described plane wave. 

2. The radius i? of a star is typically 8 orders of magnitude larger than the coherence 
length of the emitted light, and therefore coherence effects at emission time can be 
neglected: The star is a chaotic radiator. The pion source has a typical size of a few fm 
which is of the same order as the coherence length of the pion. Therefore, interference 
effects in the source have to be taken into account, as it is consistently done in the 
wavepacket model [4-6]. 

3. Photons are neutral, whereas the pion correlations (with exception of 7r°) are modified 
by the Coulomb interaction. 

The Coulomb distortion of the correlation function and the corresponding correction has 
been the topic of many recent publications [7-10], where the pions where described either as 
plane waves or as classical particles. The several existing theoretical approaches all have in 
common that they do not consider the zero point energy which emerges as a consequence of 
the finite size of the pion source. However, in a recent wavepacket calculation of identical 
and Coulomb interacting pions it was found that the correlation function does not exhibit 
any significant distortion, a fact which was explained by the strong momentum uncertainty of 
the wavepackets [0]. This argument, if valid, should also apply to the correlations of 'k~^-h~ 
pairs. Here, positive correlations were found experimentally [12|, which are interpreted as 



caused by the Coulomb attraction between both pions and equivalently used to Coulomb- 
correct the correlation functions for like-charge pions. 



It is therefore of crucial importance to verify the conclusions of Ref. [11| also for the vr"'"- 
vr" system of non- identical mesons, where the wavefunction does not have to be symmetric 
with respect to particle exchange. This system therefore displays the consequences of the 
quantum behaviour, and in particular those of the Heisenberg uncertainty principle, most 
purely, and it allows to study the modifications when the system approaches the classical 
limit. To this end, we compute the correlation function of ti~^-'k~ pairs using wavepackets for 
the pion-states. We again discuss the important role of the dispersion which determines the 
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size of the Coulomb distortion. We compare the results with experimental data and finally 
discuss the implications. 



2 The multiconfigurational procedure 

The method to solve the time dependent Schrodinger equation for 7r"^-7r~ pairs, 

H^{ri,r2,t) = ihdt'H {ri, r 2, t) , (1) 

with the Hamiltonian 

H = -^{A, + A2) + - 1 (2) 



is similar to the method used in For the initial 2-particle wavefunction we use the 

product state 

^'(ri,r2,t = 0) = V'i(ri,t = 0)-V2(r2,i = 0) , (3) 
where the single particle states are Gaussians, i. e. 

^P,ir,,t = 0) = i2nalr'/' exp (^^P, • r, - (4) 

with j G {1,2}. Rj is the initial position of the j'th wavepacket-centre, Pj the initial 
momentum and o"o the initial wavepacket- width. Transforming into center of mass (com) and 
relative (rel) coordinates, the state is rewritten as 

^'(rc, r, i = 0) = ^'com(rc) • *rei(r) (5) 

with Tc = (ri + r2)/2 and r = ri — r2, 

^com(rc) = (nal)-'/' exp (^^Pe • r, - , (6) 

and 

*,ei(r) = (4vra„2)-3/4 exp (^^P • r - ^^^^ ■ (7) 

Here, Rc = (Ri+R2)/2, R = R1-R2, Pc = P1 + P2 and P = (Pi-P2)/2. Since there are 
no external forces, the time dependent solution for ^'com is the solution for a free wavepacket, 
i. e. 

*com(rc, t) = {ns.itfr^/' exp {P, • - E,t) - • (8) 

\n 2sc{t)ao J 

Here, = P^/{2M) with M = 2m and the pion mass m = 140 MeV, finally 
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In order to evaluate the time dependent solution for ^'i.ci(r, t), an expansion into a set of basis 
functions is required. As basic components for our basis we use the same time dependent 
wavefunctions as in 111, i.e. 



y^A^x' exp U {P{t) ■ r - 9i,m,n{t)) " ^^^^^^ 



with the phase 



Qlmn{t)= / -E(t) dr + ?i (/ + m + n + 3/2) arctan 



ht 

2^ 



(10) 



(11) 



and = (r^; — R^)', and similar for the y- and z coordinates. Further, P = (Pi — P2)/2, 
is the reduced mass, /, m and n are positive integers. 



s{t) 



0"o 



1 + 



iht 



and the norm is 



Ml 



Lm,n 



(W)-3/2(V2a)-2('+-+") 



fl2) 



(13) 



(2/-l)!!(2m-l)!!(2n-l)!! 

with a = \s\ and j\\ = j {j — 2) (j — 4) • • •. We may call a state with k = l + m + n = 
"s-state", similar "p-state" for k = 1, "d-state" for k = 2 and "f-state" for A: = 3. The basis 
functions il^i,m,n{^:t) are not orthogonal for different sets {l,m,n}, but linear combinations 
of them can be found which are orthogonal [13|. We use the following orthogonal basis: 



■0000 

V'lOO, 0010, 0001 

0'2OO - "0020, 0^011, 2-0002 " 0200 " 0'O2O, 0^101, 0110 

0300 - 30120, 0201 - 0021, 4-0102 " 0300 " 0120, 0030 " 30210, 

20003 - 30201 - 30021, 40012 - 0210 - 0030, 0111- 



The expansion of ^'rei i^ow reads 



^'rci(r,t) = 5Icfc(t)04(r,t) 



(14) 



(15) 



with the c-valued coefficients 



Ck{t) = (0.^-;(r,i)|M/rei(r,t)) (16) 

and the index- vectors 1^ = {Ik^rrik^nk)- The equation of motion for the expansion coefficients 
is 

inc= (h -b^ c (17) 

with the time dependent coefficient vector c, the Hamilton matrix Hij = (0^-* |lf|0^-' ), and the 
time evolution matrix Dij = (0r|i?ii9t|0r ). 



3 



Equation (|17| ) can be regarded as the quantum mechanical part of the time evolution, 
which is necessary to account for the deformation of the wavepacket. A second, classical 
contribution is given by the motion of the wavepacket-centre. The corresponding equation of 
motion is defined as 

P = -Vn{^relir)\Vir)\%,i{r)) (18) 



with V{r) = —e /\r\. The evaluation of the matrix-elements is described in |13]. The 
numerical solution of the equation of motion was done the same way as described in |11], i. e. 
a second order integrator was used for the integration, the simulation was stopped after the 
Coulomb-potential has dropped to 2% of its original value, and the final state was Fourier- 
transformed into momentum space. In order to evaluate the 2-particle correlation function 
of the relative momentum q = |pi — P2I, 

a Monte-Carlo procedure Q was used to sample the 2-particle probability density 'P2(Pi> P2)i 
and event-mixing was used to obtain the denominator of Eq. (|T9|). As described in [|ll|, it 
is possible to generate a large number A'^, of pairs from one single momentum state. After a 
number of Ne independent events are calculated, the total number of pairs equals Np = Ne-Ng. 

3 Results 

As source function, i. e. the initial phase space distribution of the wavepackets, we applied 
the following simple parametrization: 

p{K) = {nRl)~^/'ew(=^^ , (20) 

/(P) = (2vrmr)-3/2 exp (^^^ , (21) 

i. e. the source was a static Gaussian with rms- radius y/3/2Rs and temperature T. The initial 
width (To of the wavepackets is to some extent arbitrary, but there exist certain upper and 
lower limits. Since any localization of a quantum state is accompanied with the corresponding 
zero point energy 

the observed pion energies [^4], |l^ do not allow the localization to be well below ~ 1 fm, 
which already yields -E'o(co = Ifm) ~ 100 MeV. On the other hand, the correlation radii 
which are observed in correlations of identical pions also depend on cTq, because the effective 
radius of the source is given in 1. order by 



i?eff = ^Rl + 2al . (23) 

Since the measured radii have a typical size iicorr ~ 6 fm the upper limit for do may not 
significantly exceed ~ 4 fm. 
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In Fig. I the correlation functions are shown for 3 different initial wavepacket widths. 
In one of the cases (<To = 4 fm, Rs = fm) the width of the wavepacket determines the 
source size. In a second case, ctq = 1.8 fm and i?^ = 5 fm yield the same effective source size 
ReS « 5.7 fm as in the first case. The simulation with o"o = 20 fm corresponds to the effective 
source size R^s = 28 fm, it is therefore not compatible with the measured radii in like-charge 
pion interferometry and should be considered as a pure exercise to study the dependence of 
the Coulomb effect on ctq- In all cases the source temperature was set to T = 50 MeV. The 
states were evaluated using the spd-hasis, in test simulations it was checked that the increase 
to the larger spdf -basis did not change the results in any observable way. The following facts 
are important: 

1. The enhancement in the correlation function for small relative momenta q is much 
weaker than obtained in the classical trajectory calculation (solid curve). 

2. The enhancement becomes more pronounced with an increasing value for the initial 
width (To of the states. 



These results confirm the observations made in simulations of identical pions [11|. There, 
the factual disappearance of the Coulomb signal was explained by the dominance of the zero 
point energy over the Coulomb energy. The maximum Coulomb energy in a system of 2 
Gaussian states is reached when both completely overlap, yielding (without interference, i. e. 
for nonidentical particles) 

i^S^f = . (24) 



We define the ratio 

f^^=v^;?-^'""° '''' 

as the relative Coulomb strength of the system. We obtain ^ = 1.4 • 10~^ for the cTq = 1.8 
fm case, ^ = 3.1 • 10^^ for Uo = 4 fm and ^ = 0.16 for Uo = 20 fm. Due to the smallness 
of ^, in all cases the system is dominated by the quantum dispersion. As a crosscheck, the 
fundamental constant H was reduced to 1/10 of its real value. This implies a reduction of Eq 
by a factor of 100 and, for ctq = 1.8 fm, a relative Coulomb strength ^ = 1.4, which is now in 
the classical domain. The simulation results (with source size Rs = 5 fm, stars) are already in 
good agreement with the corresponding classical trajectory-calculation for pointlike particles 
{Rs = 5.7 fm, solid curve). 

Similary, an increase in particle mass and/or charge can shift the ratio ^ towards the 
classical domain. For example, protons with o"o = 4 fm correspond to ^ = 0.2, a-particles 
with do = 4 fm correspond to ^ = 3.4. In simulations, a strong Coulomb signal was observed 
for these particles. Nevertheless, for small relative momenta, a deviation from the classical 
results remains, since the momentum uncertainty dp = Ti/{2ao) = 25 MeV/c does not depend 
on the particle mass, but reduces the momentum resolution which becomes significant for the 
rapid rise of the Coulomb signal when q approaches zero. 

4 Discussion 



In proton- |1£] or nucleus- |12] induced reactions with energies above 10 AGeV, rather strong 



positive correlations in 7r"^-7r pairs have been observed. These are interpreted as a Coulomb 
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effect — a conclusion which is supported by classical trajectory and/or plane-wave scattering 



calculations. In this work and in |11| we have tried to argue that both assumptions are 
not applicable for particles which are emitted from a source, the size of which is 2 orders of 
magnitude smaller than the Bohr-radius of 7r"'"~7r~ pairs. Because of Heisenberg's uncertainty 
principle a spatial uncertainty of size Gq transforms into a momentum uncertainty which 
allows the mutual particle interactions to become discernible in momentum correlations only 
when the parameter ^ in Eq. (^) is of order 1. The value of ^ can increase by either increasing 
the interaction strength, the delocalization of the states or the reduced system mass. For 
the TT~^-'K~ system, the corresponding parameters under normal conditions are too small. 
Therefore, it is highly desirable to find an interpretation of the ii~^-'k~ signal. Obviously, 
there are two ways out of the dilemma: 

1. The observed correlations are of Coulomb origin, and there is a yet unknown inconsis- 
tency in the wavepacket model and/or its interpretation. For example, if the pion-states 
could be initiated with an arbitrary large wavepacket width ctq, then the ratio ^ in Eq. 
( p5| ) easily could reach values above 1 and hence the classical (plane-wave) assumption 
would be justified. 

2. The observed correlations have another origin, for example strong interactions, or resid- 
ual dynamical correlations which could not be eliminated in the event mixing procedure 
during data analysis. 

With respect to the first point, i. e. the wavepacket size, one has to bear in mind that the 
pion wavefunction is interpreted as the probability amplitude to detect the pion at position 
r in case of a measurement, and there is no physical reason to assume that the pion can 
be found far out of the source already at emission time. Additionally, the short lifetime of 
the A(1232)-resonance, the scattering and reabsorption by the surrounding nuclear matter 
and, particulary, the experimentally observed small values for the effective source sizes do 
not suggest the description of the pion state by a widely extented coherent wavefunction. 
Consequently, ctq cannot exceed the size of the source, especially since that size is to be 
determined in the interferometry of like-charge pions. In order to measure an object with 
a given spatial extent Rs-, one unavoidably has to cope with the corresponding momentum 
uncertainty. 

But there exists a third and to some extent speculative possibility, which has also been 



discussed in |16|: With increasing energies, more and more long-living resonances are formed, 
which do not decay within a small volume. Pions which are produced under such circum- 
stances do not necessarily fall under the above considerations, i. e. they may approximately be 
described as plane waves. These pions could contribute to the observed 'k~^-tt~ peak, but this 
process would not work for it~^-tt~^ and 7r~-7r~ pairs emitted from the short lived and smaller 
interaction zone formed in heavy ion collisions. Consequently, the measured correlation func- 
tion of like-charge pions should be corrected by a Coulomb factor which is weighted with the 
fraction of pions produced by long-living resonances. This hypothesis does not explain why 
the observed tt^-tt^ correlations are of the observed strength, since evidently not all pions 
could have been produced far outside the interaction zone — otherwise, there wouldn't be 
any visible Bose-Einstein signal in the like-charge data. It seems that low energy data, which 
exclude the influence of long-living resonances, could help to solve the ■k~^-'k~ puzzle. 
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■ a^ = 1.8fin 
• c^ = 4fm 

A <5^ = 20fm 

^ <5^ = 1.8fm, h-^h/10 

— classical particles 
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Figure 1: The correlation function for 7r"'"-7r~ pairs with different initial localizations Gq- 
Systems with (Jo < 4 fm (solid symbols) are compatible with the results from like-charge pion 

interferomctry. The simulation with reduced Planck-constant (stars) displays the approach to 
the classical limit (solid curve). Number of generated pairs: 5- 10^ (squares), 3- 10^ (circles), 
5 • 10^ (triangles), 5 • 10^ (stars). 
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